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Abstract - We analyze the stability of a low-order coupled sea ice and climate model and extract 
the essential physics governing the time scales of response as a function of greenhouse gas forcing. 
Under present climate conditions the stability is controlled by longwave radiation driven heat 
conduction. However, as greenhouse gas forcing increases and the ice cover decays, the destabi- 
lizing influence of ice-albedo feedback acts on equal footing with longwave stabilization. Both are 
seasonally out of phase and as the system warms towards a seasonal ice state these effects, which 
underlie the bifurcations between climate states, combine exhibiting a "slowing down" to extend 
the intrinsic relaxation time scale from ~ 2 yr to 5 yr. 
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1^ Introduction. — Earth's climate is controlled by solar 
L_Z^ncrgy input, longwave radiative output and the internal 
redistribution of energy by the atmosphere and the ocean. 
Of the 5 PW excess energy input at low latitudes about 
two thirds is transferred poleward by the atmosphere and 
the remainder by the ocean. One can ascribe an energy 
OOtiudget to each subsystem and among these the polar re- 
^ions are of particular interest due to their strong influence 
^^on the planetary albedo. The ice-covered polar oceans 
dre understood to be a sensitive diagnostic subsystem for 
changes in climate, because the contact between the atmo- 
[ Sphere and the ocean is mediated by a thin veneer of sea 
J> ice that has an albedo about a factor of three larger than 
\rl that of the ocean. The ice-albedo feedback is always posi- 
tive amplifying perturbations in ice coverage and tending 
Cdto drive large scale glaciation or deglaciation [1]. Indeed, 
the magnitude of the seasonal areal change of sea ice is 
10-20 M km^ in the Arctic and Antarctic respectively, and 
hence constitutes one of the largest variations in surface 
albedo on seasonal time scales as does any component of 
the climate system. While basic physics and paleoclimate 
data show that large increases in greenhouse forcing lead 
to the decay or vanishing of the ice cover, the state of 
the art numerical global climate models have not made 
reliable projections of the substantial rate of decay of the 
minimum annual ice coverage observed by satellites [2] . 

The following facts motivate this study, [i] There are 
only about thirty years of satellite observations, during 
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the last decade of which the ice cover has decayed rapidly 
[5]. {ii) Although we understand the leading order physics 
governing the system, the state of the art global climate 
models do not reliably capture these observations. {Hi) 
We cannot use measurements to close the energy budget 
with the resolution needed to account for the observed ice 
loss P]- Thus, we are emboldened to appeal to a quantita- 
tive stability analysis of a qualitative, but observationally 
consistent, theoretical treatment to assess the nature and 
the characteristic times scales of the transitions we may 
expect. 

Recently Eisenman and Wettlaufer [3] (here referred to 
as EW09) developed a low-order theory that describes the 
seasonal evolution of the energy state in the Arctic Ocean. 
The theory couples the poleward atmospheric heat trans- 
port to the surface radiative heat balance through a gray 
body atmosphere and to the base of the ice via ocean 
heat flux. Over climatological time scales the evolution 
of the ice state is solved as a two-moving boundary prob- 
lem; perennial ice grows or thins seasonally according to 
the energy flux balances across its thickness. As the green- 
house forcing increases ice is lost in summer but returns in 
winter. Finally, under further warming this seasonal state 
gives way to a perennially ice free state via a saddle-node 
bifurcation. In analogy with a first order phase transi- 
tion, it is basic to such a bifurcation that a hysteresis 
exists between the two perennial states; in this case ice 
covered and ice free. Thus, the perennial ice cover only re- 



p-1 



W. Moon and J.S. Wettlaufer 



turns after the chmate cools to a second bifurcation point 
where the perennially ice-free Arctic would be sufficiently 
cold — colder than the original climate state at which the 
ice disappeared — to freeze. Here, we analyze the stabil- 
ity of the seasonal steady states of this model to extract 
the principal stabilizing and destabilizing effects and their 
associated response time scales. 

Response Theory of the Dynamical Model. 

Here, we summarize the theory of EW09 and refer the 
reader to their Appendix for the full derivation. The state 
variable E is the energy (with units W m~^ yr) stored in 
sea ice as latent heat when the ocean is ice-covered or in 
the ocean mixed layer as sensible heat when the ocean is 
ice-free, viz.. 



E 



—Lihi E < [sea ice] 

CmiHmiTmi E >0 [occau] 



(1) 



where Li the sea ice latent heat of fusion, hi its thick- 
ness, Cmi is the specific heat capacity of the ocean mixed 
layer, Hmi is its depth and T^i its temperature. Ignoring 
colligative effects, the temperature T{t, E), determined by 
energy balance across the layer, is measured relative to the 
freezing point Tm as 



T{t,E) = -n 



hUlE - FT{t) 



(2) 



where the ramp function is TZ{x > 0) = a; and TZ{x < 
0) = 0, the thermal conductivity of the ice is ki, and the 
radiative quantities Fo and Exit) are discussed below. 

The Beer-Lambert law of exponential attenuation of ra- 
diative intensity with depth in a medium requires a treat- 
ment of the dependence of the surface albedo with E. 
While this is done with parsimony it is physically realistic; 
one uses a characteristic ice thickness ha for the extinction 
of shortwave radiation as follows 



ami + Oti , ami 



-tanh 



E 

Ljhn 
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which controls the fraction, \ — a{E), of the incident short- 
wave radiation Fs{t) absorbed by the ice [3]. 

The evolution of the state of the ice (or ocean) cover is 
determined by the balance of radiative and sensible heat 
fluxes at the upper surface, Fd — FT(t)T(i, i?), the upward 
heat flux from the ocean Fb^ and the fraction of ice ex- 
ported from the domain vqTZ{~E) ('^ 10% yr^^) through 
a first order nonautonomous energy balance model as 



with 



f[t,E) =Fd- FT{t)T{t,E) + FB+vaTZi-E), 
where 

FD{t,E) = [1 - aiE)] Fs{t) - Fo{t) + AFo. 



(5) 



(6) 



The term Ed - FT{t)T[t,E) is thought of as the differ- 
ence between the incoming shortwave radiation at the 
surface [1 — a{E)] Fs{t) and the outgoing longwave radi- 
ation (oc T*), augmented here by sensible and latent heat 
fluxes and an additional amount associated with green- 
house gas forcing AFq. The Stcfan-Boltzmann equation 
for outgoing longwave radiation is linearized in the devi- 
ation of the surface temperature from the freezing point 
as Fo{t) + FT{t)T{t, E) which is known to be a reason- 
able approximation [3l|4]. Finally, the seasonally varying 
values of Fo(t) and Exit) are determined from an atmo- 
spheric model incorporating observations of Arctic cloudi- 
ness, atmospheric transport from lower latitudes and the 
meridional temperature gradient [3|. 

The seasonal steady states (Es), or periodic points, of 
Eq. (m were described heuristically in the introduction 
and examples arc shown Fig. 3 of EW09, wherein sensi- 
tivities to greenhouse forcing and other parameters were 
studied. Here we assess the nature of the stability of these 
periodic points and their relaxation time scales. 

Consider E{t) = Es{t) + ^(i), where \^{t)\ < \Esit)\ 
and thus 
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dEs 
dt 



d^ 
dt 



= fit,Es+0-fit,Es 
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t (7) 



E=Es 



Now, because dEs/dt = f{t,Es) we have 

d^^dl 
dt dE 



c ^ am, 



E=Es 

which has an exact solution 



^(T) = ^(O)exp 



a{s) ds 



me'' 



(8) 



(9) 



written here for the annual cycle of the system T. There- 
fore, an unstable (stable) periodic point of Eq. Q has 
7 > (7 < 0) with a perturbation relaxation time scale 
of divergence from (convergence to) a particular periodic 
point determined by the competing effects embodied in 
a{t), itself depending on whether there is ice, E < 0, 
Foit) < 0, or ocean, E > 0. li Foit) > the ice ab- 
lates downward from the surface while T{t, E) is pinned 
at zero. Accordingly, we examine the structure of a{t) in 
the two thermodynamic regimes; the ice covered and ice 
free states. 

Ice Covered State: E < 0, F£i{t) < 0. Generally we 



find 
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^ aiAjt) 
~ Fs{t) ' 
(11) 
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dT{t,E) 
dE 



E=Es 



and 

aiAjt) 
" -hL,/Es + FTit) 

[-kiLjEs+FritW Es^ 
^ aARjt) + aLwjt) 
Frit) 

Therefore, we write Eq. (fTU|) as 

a{t) = aiA{t) + aAR{t) + aLw{t) + aEx{t), 



(12) 



(13) 



where ajAit) describes the ice-albedo feedback, aAR{t) 
and aLw{t) the albedo and longwave responses, and 
dExit) = —vq the ice export, each of which we describe 
below. Note that the derivative of T{t, E) with respect 
to E in Eq. ([T2|) depends on the sign of Foit). When 
FD{t) > and the ice is ablating dT{t,E)/dE \e=Es= 
so that a(t) = ajAit) + aExit). 

Ice-Free State: E > 0, Foit) > 0. The response 
of the system is limited by the radiative balance over the 
ocean mixed layer and hence we have 



ait) 



Frit) 



(14) 



Dissecting the Stability Parameters; a{t) &: 7 . 
First we describe the climatological evolution of the direct 
response rate a(i), for it underlies the seasonal evolution 
of the ice state. Next in more detail we analyze the in- 
tegrated influence of those dynamics on the perturbations 
to the system viz., 7. Note that the dominant underly- 
ing physics influencing stability is laid bare by this simple 
analysis, but is included in more complex numerical mod- 
els (see e.g., [5] and refs therein). 

The evolution of the response rate a{t) and its com- 
ponents as described by Eq. ((T3| are shown in Figure [1] 
along with the seasonal steady state solution Es{t) (blue 
lower curve) described by Eq. and Foit) (green lower 
curve) given by Eq. The net radiative forcing and 

the surface temperature determine which components of 
a{t) are dominant. During winter, when FD{t) is nega- 
tive, ULwit) dominates, whereas during summer qja dom- 
inates. As greenhouse forcing AFq increases the dynamics 
of a{t) change and we see that when AFq = 19.0 W m~^, 
although the steady state solution still represents a peren- 
nial ice state, during summer sea ice becomes sufficiently 
thin that sea ice albedo feedback nearly drives an insta- 
bility in the ice cover. Indeed, the principal difference 
relative to the AFq = 10.0 W m~^ case is the magnitude 
of ajAit) during summer which controls a{t); starting in 
April the ice starts to thin and by May its thickness lies 
in the range where the ice albedo feedback controls Es (t) . 
However, by September Fult) turns negative and thin ice 
grows rapidly. By AFq = 21 W m~^ the seasonal ice state 
appears. Whilst these dynamics provide a detailed picture 
of the seasonal variation of the stability of the system and 
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Fig. 1: The response rate a{t) of Eq. H13|) and all of the contri- 
butions defined therein as a function of month for three values 
of of greenhouse gas forcing AFo of 10 (top), 19 (middle) and 21 
(bottom) W m^^, along with the evolution of the energy Es{t) 
(lower blue curves) described by Eq. and the radiative forc- 
ing Foit) (lower green curves) as given by Eq. ((G]). The solid 
black curve is the net response rate a{t) and the individual con- 
tributions are denoted by the subscripts of Eq. (|13p . We see 
the strong seasonal dependence of the dominant contributions 
to the stability; destabilizing ice-albedo feedback (a/A > 0) and 
the stabilizing influence of heat conduction driven by radiative 
loss at the ice surface {aLw < ). 



their influence on the seasonal steady state solutions i?s [t) 
as greenhouse forcing increases, it is the seasonally aver- 
aged contributions embodied in 7 that most succinctly 
describe the overall stability. 

Following the same notation as in Eq. (fT5)) Figure [5] 
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shows 7 and its constituents as a function of the green- 
house forcing AFq. The range of AFq shown spans the 
perennial and seasonal ice states of the system. In the fol- 
lowing paragraph we summarize the principal results, and 
then discuss the contributions to 7 in turn. 
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Fig. 2: The response of the ice cover 7 to perturbations as a 
function of greenhouse gas forcing AFq. When 7 < (7 > 0) 
the periodic points Es of Eq. are stable (unstable), and 
their relaxation time scale is |1/7|. The solid black curve is 
the entire response which for small AFq is dominated by the 
longwave contribution {yLW < 0; blue dashed curve). As the 
ice thins the ice albedo feedback (7/A > 0) figures more promi- 
nently (red dashed curve) and the stability is dominated by the 
competition between the destabilizing ice-albedo feedback and 
the stabilizing influence of heat conduction driven by radiative 
loss at the ice surface. The local maximum at AFq — 19.8 
W m~^ is the transition from the perennial to the seasonal 
(winter only) ice state and the second increase at AFq = 22.2 
W m~^ denotes the saddle-node bifurcation to a seasonally ice 
free state. The small differences between these thresholds and 
EW09 are because we use a{E) rather than a; in Eq. Q. 

For greenhouse forcing up to AFq ^ 16 W m"^, the 
perennial ice is sufficiently thick that the ice-albedo feed- 
back plays a minor role and longwave radiative control 
of ice growth in winter dominates the stability with a 
rather rapid relaxation time of a few years. As the cli- 
mate warms and the ice thins the stabilizing influence of 
heat conduction acts on a much more rapid relaxation time 
scale (thin ice grows more rapidly than thick ice [5]) but 
so too docs the destabilizing ice-albedo feedback governed 
by Eq. ((TT|) . These competing effects lead to an overall 
increase in the relaxation time with greenhouse forcing of 
^ 5 yr until AFq ~ 19.8 W m~^, ice is lost in summer, 
and reentrant longwave control of ice growth in winter 
dominates the rapid stability of this state of the system. 
Finally, when AFq > 22.2 W m~^, the ice is too thin to 



thwart the ice-albedo feedback, it does not recover in the 
polar night, and is thus lost entirely. 

The first contribution, aiA{t), describes the fluctuation 
in the surface shortwave radiative flux associated with a 
perturbation in ice thickness. The fundamental operation 
of the ice-albedo feedback is laid bare by the fact that 
aiA{t) is always positive; amplifying all perturbations in 
ice thickness. Eq. [TT]is rather transparent in the display 
of the albedo contrast as a function of the ice thickness. 
Although Fs{t) — during winter, the annual manifesta- 
tion of the ice-albedo feedback is jia = Jq aiA{s)ds > 
and, as AFq increases, this drives the transition to both 
the seasonal ice and ice-free states. Such behavior is due 
to the fact that jia only begins to play a controlling role 
in the state of the system when h < ha. Importantly, a 
value for of 0.5 m is realistic given the optical proper- 
ties of sea ice. However, an unphysical value of even half 
of this results in unrealistic transitions in the state of the 
system; albedo is a material property. 

The Stefan-Boltzmann law underlies the principal sta- 
bilizing influence, understood from rudimentary-radiative 
steady state-studies of the greenhouse effect. Increasing 
the flux Fo{t) increases the surface temperature T{t,E) 
and hence the outgoing longwave flux via the Stefan- 
Boltzmann law. In this coupled model it is prudent to 
take care in interpreting this feedback. Here, any change 
in T{t, E) has an immediate influence on the heat con- 
duction through the ice but the atmospheric meridional 
heat flux is proportional to the difference between the ice 
surface temperature and that at the lower latitudes. For 
example, when Fs{t) increases during early summer, both 
T{t, E) and the associated outgoing longwave radiance in- 
crease, but the atmospheric meridional heat flux decreases. 
Thus, to leading order these two effects have the same sign. 
Whereas, in winter Fs(t) = and any anomalously thin 
ice grows more rapidly due to the basic tenets of heat con- 
duction driven by longwave loss at the surface. Combined 
these processes arc captured by aAR{t) and a^wit). It is 
evident that because aAR^t) = aiA{t)/[—kiLi/ Es+Fxit)], 
this term describes the stabilizing effect of heat conduction 
and outgoing longwave radiance weighed against the per- 
turbation in incoming shortwave radiance associated with 
the fluctuation in the surface albedo. Whence, we call 
this "albedo response". Whilst jab. < 0, its magnitude 
is small because T{t,E) = and hence aAR{t) = when 
the ice is ablating during summer, and aAR{t) = dur- 
ing winter because Fs{t) = 0. Therefore, the stabilizing 
action of aAR{t) is conflned to the spring and fall. 

The principal influence embodied in a^w is the growth 
of ice by longwave cooling during winter; perturbations 
associated with ice-albedo feedback when Fs{t) > are 
strongly compensated for during the polar night. The 
thinner the ice at the end of the summer the greater, and 
more rapid, the stabilizing growth response. While the 
ice ablates in summer T{t,E) = and hence ULwit) van- 
ishes. However, it describes the dominant stabilizing effect 
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in the system insuring that ^lw < with |7lw| increasing 
as AFq increases and the ice thins. 

Finally, the rate of sea ice export out of the Arctic 
Ocean, controlled by wind and ocean forcing, is repre- 
sented by "fEX, and when the ice vanishes the longwave 
response is shown as jni- 

The dominant competition determining the stability of 
the periodic points of the theory, and the associated re- 
laxation rate as AFq increases, is between the ice-albedo 
feedback aiA{t) and the nonlinear response of ice growth 
associated with longwave loss in winter aLw(t)- These two 
effects act in opposite seasonality and hence provide an in- 
tuitive picture for how one state can emerge from another 
as the persistence of the ice cover is enhanced/suppressed 
due to fluctuations in forcing. Moreover, once AFq has 
increased beyond about 19.8 W m~^ this competition 
demonstrates the origin of the robustness of winter ice. 
Minor manipulation affords some insight into the main 
stabilizing effect viz.. 



a-Lwity 



FTit)f Es 
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hsFrit) 
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where T is the perturbation in surface temperature and 
is the steady state ice thickness. This form of aLw{t) dis- 



plays the stabilizing effects 



and M^; the longwave 



heat loss controlling heat conduction, and the fraction of 
radiative forcing stored as latent heat in the ice respec- 
tively. The former effect is basic to heat conduction; for 
a given longwave radiative forcing thinner ice grows more 
rapidly than thicker ice. Thus, as the climate warms and 
AFo increases, hs decreases leading to a faster response 
for the stabilization of perturbations. Considering annual 
mean values (see Table 1 of EW09) for demonstrative esti- 
mates f'^^*2"^ sa 4.2 for 3 m thick ice under no greenhouse 

kil/hs 

forcing (AFo=0) and about 1.4 for 1 m thick ice with AFq 
= 19.8 W m~^ (where perennial ice is lost). Inspection of 
Eq. (|n]) shows that the shortwave absorption at the sur- 
face decreases as the ice thickens, rapidly saturating for 
h > ha and, for a given ice thickness, increases with AFq. 
For the same conditions given above, decreases from 
about 0.5 to 0.3 yr with increasing AFq. Taken together 
this demonstrates the contributions to the more rapid sta- 
bilizing growth response of thin ice in a warming climate 
as |aLvi/(i)~^| decreases from ~ 3.5 to 1.2 yr. 

Conclusions. — As seen in figure [2] the overall stabil- 
ity and the net response rate of the periodic points of Eq. 
Q to perturbations depends on all the contributions to 7, 
but it is dominated by the competition between the desta- 
bilizing ice-albedo feedback jja > and the stabilizing in- 
fluence of heat conduction driven by wintertime radiative 
loss at the ice surface jlw < 0. The order of magnitude 
estimates given above for the latter effect demonstrate the 
general stabilizing behavior that thin ice responds much 
more rapidly than thick ice to perturbations in radiative 
forcing. Because ice with h < ha absorbs more shortwave 



radiation than does thick ice (Eq. [3]), the destabilizing ice- 
albedo feedback also operates much more rapidly as AFq 
increases and the ice thins. Thus as AFq increases while 
these principal, seasonally out of phase, contributions to 
the competition which governs the stability of the system 
(7/yi and 7lh/) each operate on shorter time scales, their 
difference decreases thereby producing a collective state 
that responds more slowly 5 yr as AFq — > 19.8 W 
m~^). Whereas under conditions of small AFq more rep- 
resentative of the present climate, we find time scales 7^^ 
of 2 yr. Hence, a recent study of the simulated recov- 
ery using an atmosphere-ocean general circulation model 
is relevant in this regard. Tietsche et al., (2011) numer- 
ically prescribed ice-free summer states at various times 
during the projection of 21*** century climates and found 
that ice extent typically recovered within several years. 
An important implication of their results is that while the 
ice-albedo feedback may drive the transition to an ice free 
summer, the stabilization and recovery of the system can 
be rapid. Not only arc the recovery processes described in 
this compact theory ostensibly the same as in the complex 
treatment in Tietsche et al., [5] but the stability analysis 
presented here shows that the time scales are commen- 
surate. Moreover, the underlying processes and their in- 
trinsic time scales are immediately accessible within our 
framework. 

We close by noting that such a stability analysis lays 
the foundation for a rigorous study of the role of stochas- 
tic forcing in the climate state of the system [T] . This is 
because while the transitions between the periodic points 
of Eq. (|4]) are essential to understanding the basic physics 
of the system [3], the nature of these transitions will de- 
pend intimately on the response time scales (7"^) to fluc- 
tuations. Because of the secular trends in 7/^ and ^lw as 
the climate warms, the detailed interplay between these 
competing effects in the presence of noise comprises a sep- 
arate detailed study. 
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